Single-hole dynamics in the t-J model 
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The quasi-particle weight of a single hole in an antiferromagnetic background is studied in the semiclassical 
approximation. We start from the t-J model, generalize it to arbitrary spin S by employing an appropriate 
coherent state representation for the hole, and derive an effective action for the dynamics in the long-wavelength 
low-energy limit. In the same limit, we find an expression for the single-hole Green's function which we 
evaluate in an 1 / S expansion. Our approach has the advantage of being applicable in one and in two dimensions. 
We find two qualitatively different results in these two cases: while in one dimension our results are compatible 
with a vanishing quasi-particle weight, this weight is found to be finite in two dimensions, indicating normal 
quasi-particle behavior of the hole in this last case. 
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I. INTRODUCTION 

Although more than 30 years have elapsed since Brinkman 
and Ricetl published their pioneering work on the propagation 
of a hole in a quantum antiferromagnet, this subject still at- 
tracts a lot of attention. Central to this problem is the question 
whether the motion of a hole in a strongly correlated system 
is coherent or not, i.e., whether the hole can be considered 
as a quasi-particle in the sense of Landau's Fermi-liquid the- 
ory, or whether, as was first suggested by AndersontB, it is 
more appropriate to view the hole as a constituent particle of 
a Luttinger-like liquid. For a long time, the interest in this 
question has been fueled by the expectation that an accurate 
description of the hole dynamics in a strongly correlated elec- 
tron system would be the crucial first step towards an under- 
standing of the physics of the cuprate superconductors. More 
recently, this interest has been revived, since the dynamics of 
holes in one- (ID) and two-dimensional (2D) antiferromag- 
nets has become experimentally accessible in angle-resolved 
photoemission spectroscopy (ARPES)ffl. The key quantity 
that can in principle be obtained from the results of ARPES 
is the spectral weight Z at the ground state energy Eg of the 
systems. A nonvanishing spectral weight, Z / 0, implies 
that the hole is a quasi-particle, e.g., a spin-polaron consist- 
ing of the bare hole with a spin polarisation cloud around it. 
On the contrary, a vanishing Z signals that the hole causes 
a global rearrangement of the original ground state such that 
the exact single-hole wave function becomes orthogonal to 
the wave function of a bare hole in the antiferromagnetic 
ground state. ARPES data for the one-dimensional copper 
oxide chain compound SrCuC<2, see Ref. ^ and for the two- 
dimensional layered compound Sr2Cu02Cl2, see Ref. [I[ show 
that the hole dynamics differ significantly between the ID and 
the 2D case (see the discussion in Ref. g). The features of the 
spectra of the ID compound fit rather well with the concept of 
spin-charge separation that has emerged from BethepAnsatz 
and bosonization studies of the ID Hubbard modeltUffl and 
from the exact. solution of the t-J model at the supersym- 
metric poinoEl Theoretically, spin-charge separation im- 
plies that for long wavelengths, the Hamiltonian of a hole can 



be decomposed into two commuting parts, H = Hh + H s , 
where both the holon naEt-Ef/j and the spinon part H s , are free- 
fermion HamiltonianstlljO. As a consequence, the single- 
hole Green's function takes the form 



G(p,r) 



dQ 
2tt 



G h (p-Q, T ) Z(Q,t) 



(1) 



Here, Gh(k,r) is the propagator of a free holon, whereas 
the function Z(Q,t) is entirely determined by the spinon 
dynamics; Q is the spinon momentum. While the holon 
Green's function Gh is a free-fermion propagator indeed, 
Im Gh(k,oj) = tt5(oj — Eh{k)) with £/,(&) the holon kinetic 
energy, the spinon contribution Z(Q, r) is a highly nontrivial 
singular function. Technically, this last feature has its origin 
in the fact that the original hole operator Ci a does not sim- 
ply factorize into a holon and a spinon operator, Cj CT ^ h\si. 
Rather, in the representation of a a by h\ and &j, an additional 
phase factor cxp(i8 s ) is needed, where the phase 6 S depends 
on the spinon operators Sj-nt^. This phase factor accounts 
for the phase-string effectLiOLl: as the physical hole moves 
through the (quasi)-ordered spin background, it trails a string 
of overturned spins behind it. Using the phase-string picture, 
Suzuura and NagaosaLU have obtained a qualitatively correct 
result for the single-hole spectcaL function in ID. We note, 
however, that Sorella and ParolaBo have derived their detailed 
results for the spinon function Z(Q, r) in ID without making 
use of a phase-string picture. Thus in ID, the different ap- 
proaches lead to the same conclusion: there is spin-charge 
separation in this case, and the hole Green's function shows 
Luttinger-liquid-like behavior. A recent quantum Monte- 
Carlo (QMC) simulation of the single-hole dynamics in the 
one-dimensional t-J model has confirmed these resultstfl. 
In two dimensions, the situation is less clear. Many stud- 
ies of the hole dynamics in the two-dimensional t-J model 
have been based on an effective Hamiltonian which describes 
the hole as a spinless fermion that is coupled linearly to the 
magnon excitations of the antiferromagnetic background. In 
obtaining the hole Green's function, the coupling between the 
hole and the magnons is then treated in the self-consistent 
Born approximation (SCBA), and a normal quasi-particle 
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peak is found in the sing le-hole spectral functionEHliEEO. It 
must be emphasized that the SCBA is a perturbative method 
so that the quasi-particle picture is inherent in this method. 
However, a finite quasi-particle weight of a hole has also been 
found-in a non-perturbative large-spin study_jOf the 2D t-J 
modelEl. Moreover, a recent QMC simulatiorcJ lends support 
to these results in a wide range of the model parameter J/t 
and of the momentum k of the propagating hole. But, remark- 
ably, two theoretical predictions by SorellaEl concerning the 
behavior of the quasi-particle weight Z(Q) for J/t = 2 at the 
antiferromagnetic wave vector Q = (yr, ir) are less well con- 
firmed by the simulations in Ref. |T} Apart from this seem- 
ingly minor discrepancy between QMC simulation and cer- 
tain theoretical predictions, there is further work which casts 
doubt on the validity of the quasi-particle picture for a hole 
in the two-dimensional t-J model: results from a high tem- 
perature series expansion for the momentum distribution of 
the particles in the t-J model show a violation of Luttinger's 
theorerro. This means that the ground state of the 2D t-J 
model is not connected adiabatically to the ground state of 
the noninteracting model as would be the case for a normal 
Fermi-liquid consisting of-,ouasi-particles. More directly, it 
is claimed by Weng et a/.l!3E2tJ that the phase-string picture, 
which explains the vanishing of the quasi-particle weight in 
ID, applies in 2D too and has the same effect there. De- 
spite this apparent failure of the quasi-particle picture, these 
authors suggest in Ref. |24] that the spectral features seen in 
ARPES can be explained satisfactorily in the framework of 
the t-J model by their theory. However, in a very recent QMC 
study Mishchenko et al£5 contradict this suggestion, as they 
confirm the quasi-particle picture already found by Brunner 
et alB. 

In view of these contradictions between numerical and ana- 
lytic analyses, we propose a new analytic approach to the cal- 
culation of the quasi-particle weight Z of a single hole in the 
t-J model. Our method has the advantage of being applicable 
in arbitrary space dimensions and is thus suitable to detect dif- 
ferences between the ID and the 2D case. In ID, our results 
are compatible with a power-law decay of the weight Z as the 
system size L increases to infinity, Z ~ L~ 2X . On the con- 
trary in 2D, this weight remains finite in the thermodynamic 
limit indicating that the quasi-particle picture is valid in this 
case. 

The paper is organised as follows: Section II contains the 
derivation of the effective action of a single-hole in an antifet- 
romagnetic background with the method of coherent statestJ. 
In the later parts of the paper, we wish to employ a semiclas- 
sical expansion, i.e., an expansion in powers of the inverse 
spin length 1/S. For this reason, we generalise the conven- 
tional t-J model in which every lattice site is either occupied 
by a spin | particle or empty, to the case where the lattice 
sites are occupied either by an object with arbitrary but fixed 
spin S (particle) or by an object with spin S — h (hole)El 
We implement this modification by an appropriate choice of 
the coherent-state representation for the particles which at the 
same time satisfies explicitely the constraint of forbidden dou- 
ble occupancy. In describing the local coupling of the hole to 
the spin degrees of freedom, it turns out to be crucial to use 



local SU(2) fields as the basic variables of the effective ac- 
tion instead of working with the local sublattice magnetiza- 
tion. In Section EH, we employ the effective action of Section 
II to obtain a path-integral representation for the single-hole 
Green's function. Our aim is to extract from this function an 
expression for the quasi-particle weight Z of the hole which 
exhibits explicitly the dependence of Z on the system param- 
eters, in particular its dependence on the linear system size L. 
To achieve this goal, we have to take recourse to the semi- 
classical expansion. In Section IV, we present and discuss our 
final results. Technicalities of the developments in Sections II 
and III are deferred to Appendices A and B. 



n. PATH-INTEGRAL REPRESENTATION 

In this Section, we derive an effective action for the motion of 
a single hole in an antiferromagnet. The underlying Hamilto- 
nian is that of the t-J model, 



H = -t J2 P (ci.c r , a + h.c.) P 

<r,r> 
a 

+J J2 Sr'Sr' ■ 



(2) 



<r,r'> 



Here 
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and 



are the creation and annihilation opera- 



tors for fermion states at site r with spin projection a, and 
S r = \c\ a cr a pC r g is the spin of the fermion. P projects 
onto states where each lattice site is either empty or singly 
occupied. < r, r' > denotes nearest neighbor sites on a hy- 
percubic lattice with lattice constant a in D dimensions. The 
main difficulty in analysing the properties of (||) lies in the 
strong interactions induced by the exclusion of doubly occu- 
pied sites. We deal with this difficulty by using an appropriate 



set of coherent states (cf. Refs. 28 



which takes these con- 
This leads to the following 



straints into account explicitely 
path-integral representation of the partition function Z: 



Z = 



f£ dr {M8Ju,) + (u\H\u)} 



(3) 



The coherent states \u) are introduced and discussed in the 
Appendix^. They are parameterized by two fields; the anti- 
commuting fields rj r , rj* which describe the hole, and the com- 
muting field <7 r € SU (2) which determines the orientation of 
the spin (see Eq. (^|) below). As is detailed in the Appendix, 
our coherent states are constructed such that they allow us to 
generalise to arbitrary spin S: Each lattice site is occupied ei- 
ther by an object with spin S (particle) or an object with spin 
S— | (hole)EJ. Then, the hopping term t in Eq. (g) exchanges 
spin S objects and spin S— i objects on neighbouring lattice 
sites. The corresponding matrix elements are worked out in 
the Appendix |a|. In the parametrisation by rj and g, the term 
entering the kinetic part of the action in Eq. takes the form 
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In our SU (2) description, we introduce the 5[7(2)-field 
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FIG. 1: Partitioning of the square lattice into plaquettes 



5rJ n + '/> 



(4) 



while in the Hamiltonian part 



\n\u) =J2 { 2St 



I 1 



J 



+ -(2S - - n r • n r ^ . (5) 



In the case S = 5, everything reduces to the ordinary f-J 
model. The unit vector 



n r = ^ Tr{cr g r <J z gi} 



(6) 



which occurs in the Heisenberg term in Eq. |] points in the 
direction of the expectation value of the spin. The magni- 
tude of the spin is S, or S — h, depending on the presence 
of a hole which is expressed by the occupation number r)*ij r . 
The amplitude of the hopping term ij*r/ r , in Eq. (|j) is mod- 
ulated by the actual state of the background spin, and this 
constitutes the main coupling between the hole- and the spin- 
subsystems. Note that, while the Heisenberg part depends on 
g r only through the bilinear term (^), the SU (2)-fields g r at 
the individual lattice sites are needed to express the hopping 
term (see the related discussion of symmetries of the action in 
Ref.g§). 

Proceeding towards the derivation of an effective action, we 
now divide the lattice into plaquettes each of which contains 
a pair of neighbouring sites of the two sublattices. We label 
these with x, the lattice points of the A-sublattice. Then, the 
plaquette x contains the fields g x and g' x = g r +ae x , on the A- 
and B-site, respectively, see Fig. |l[ 

In the conventional description of a two-sublattice antiferro- 
magnet, one would introduce now the magnetisation M = 
Sa + and the sublattice magnetisation N = — of 
each plaquette as new fields and integrate subsequently over 
M. 



3x fx ; 



(7) 



we choose the SU (2) fields g x and M x as our new variables 
and, finally, we will integrate over the fields M x . A useful 
parametrization of M in terms of three real variables a, 3? u, 
and 3 u is 



M 



u 1 
i u* 



(8) 



(here, the label x has been suppressed). Then, for an arbitrary 
phase a, the point u = corresponds to the antiferromagnetic 
alignment of the spins within the plaquette, since 



n x n x = i Tr {M±a z M x a z } 



= -l + 2|u s 



+ 0(ui) 



(9) 



Next, we express all terms in the action by M x and <? x ; we 
start with the Hamiltonian part. The vector between near- 
est neighbor plaquettes is denoted by S — 2ae x , ae x ± 
ae y , ... (or 0), see Fig. [j]. Then, the hopping term between 
two plaquettes in Eq. (||) contains the expression 



(10) 



Since we are aiming at a gradient expansion of the £{7(2) 
fields in the Hamiltonian part of the action, we replace 
s in the hopping term in Eq. (B) by 



fx 9: 



x+<5 



.9 X 9 x +s =■ 1 + iS ■ A x , 
where A x becomes the spatial derivative 

A x = -i gi V g x . 

in the continuum limit. The Heisenberg term in Eq. 
also be conveniently expressed by M x and A x : 



(ID 



(12) 



n x+5 ' n x 

= g Tl " {- 9 x+<5 9 * ° z 9 * 9 *+ S CT2 } 

= X - Tr { [1 + iS ■ A x ] f M x o z M x [1+iS- A x ] u z } (13) 



In the term entering the kinetic part of the action, Eq. (Q), we 
introduce the time derivative 



91 d T 5 X 



(14) 
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Then, we have on the B-sites, 



£ 8 T g' x = M x iA T x Ml + M x d T Ml 



(15) 



Next, we collect all terms and expand the expressions in the 
exponent of the functional integral for Z, Eq. (^|), up to the 
second order in u x . Finally, we perform the Gaussian inte- 
gration over u x and identify the effective action A = jf drC 
as the exponent in the resulting path-integral. The remain- 
ing integration fields are defined on the A-sublattice: on each 
plaquette there is one commuting field g x £ SU (2) for the 
remaining spin-degrees of freedom, and there are two anti- 
commuting fields r/ x and ij x = e 4Qx r) r+a ~ on each pla- 
quette which describe the holes. Note that the phase a x 
has been absorbed in the redefinition of the hole field on 
the B-sublattice; it disappears from the effective action. We 
work in the leading order in the gradients, A x and A x so 
that the effective Lagrangian C is accurate in second or- 
der in these gradients. C can now be written as a sum of 
terms, each representing the contribution of a fixed number 
of holes. We retain only the zero- and one-hole contributions: 

C = -2MDJS 2 + £ nlam + £i_hole H ■ The first term 

is the energy of the classical antiferromagnetic state of a sys- 
tem with Af plaquettes. The fluctuations of the spins in the 
absence of holes are governed in the continuum limit by the 
well-known non-linear er-model for the unit vector n x 



C 



nlcrm 



2a 



D 



4DJ 



il + JS 2 a 2 J2 



S f d u x . 



(16) 



The first two terms in Eq. ( |T6[ ) are temporal and spatial fluc- 
tuations, respectively (h x = d T n x ). The last term is the 
famous topological term which yields a zero mass for the fluc- 
tuations for half-integer spins in one dimension; for D > 1, 
this term is not effective (see, e.g., the discussion in Ref. |30| ). 
We wish to emphasise that we obtain here exactly the same 
result, Eq. (|l6|), as in the conventional approach in which one 
uses the magnetisation M = Sa + Sb and the sublattice mag- 
netisation as N = Sa — Ss as plaquette variables and then 
integrates over M. However, our approach which uses the 
SU (2)-fields g x as basic variables turns out to be superior 
when the steps leading to C n i am have to be performed for the 
hopping term of Eq. (|5(). From the procedure described above, 
we get 



A-hole — 



[?7 x +<5 ^ x +(5, x Vx + ^ x Q x +(5 Vx+S 

x,<5 



(17) 



where h a .^ are the matrix elements of the effective Hamilto- 
nian for the hole in the background of the g x -fields. The diag- 
onal elements are expressed in terms of the vector n x (bilinear 
in g x and g x , see Eq. (n)) as follows 



^^nl - ^(n x+5 - n x+<5 ,) 2 



8DJS x 16D 



8,8' 



8DJS x 16D 



8,8' 



-i i n x • ri x x i ^(n x+5 - n x ) . (18) 



The off-diagonal elements describe hopping from the A- 
sublattice to the B-sublattice and vice versa. They can only 
be expressed with the aid of the SU (2)-field g x itself: 



,0 b 



h ba 



x , x +(5 



-it 
~DJ 

it 
~DJ 



gi d T g x + JSJ2 (di+s - al+8') 9, 

8' 

9* d T g x + JS ^ .4 (g x+s - g x+s ,) 



I [ 

Tl ' 
(19) 



The preceding equations complete the derivation of the effec- 
tive action for a hole in the background of the spin-fields g x . 
The approximations made so far were (i) the Gaussian integra- 
tion over the variables u x which describe the deviation from 
the antiferromagnetic order within a plaquette and (ii) the gra- 
dient expansion in temporal and spatial derivatives. Thus, the 
Euler angles defining the fields g x should be smooth in space 
and time, but there is no limitation on their variation over large 
distances in space or time. 

We remark that our effective one-hole Lagrangian, Eq. (|T7|), 
differs significantly from the effective Lagrangian that has 
been derived by Shankar in Ref. |3l| , [32| for the same physical 
situation. Partially, the difference stems from the fact that the 
work in Ref. ^l],|2]is based on a generalised t-J model which 
includes besides nearest neighbour (intersublattice) hopping 
also next-nearest neighbour (intrasublattice) hopping. Fur- 
thermore, however, terms that describe intersublattice hop- 
ping, i.e., the terms in the last line of Eq. ( 17 ), are negelected 
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in Ref. p 1|J32| . It is argued that intersublattice hopping pro- 
cesses are necessarily accompanied by spin fluctuations and 
are therefore supressed in a situation with strong short-range 
antiferromagnetic order. In our case of a t-J model with 
purely nearest neighbour hopping, the neglect of the corre- 
sponding terms in the effective Lagrangian would leave us 
with a static hole which is clearly not an appropriate descrip- 
tion of the physical situation which we intend to consider. 



HI. GREEN'S FUNCTION 

In order to obtain the quasi-particle weight Z of a hole in an 
antiferromagnetic background, we study the Green's function 
of a single hole in the t-J model 



arrive at the formal result 



G^(r-r';r) = ( 4,» c r , CT ,(0) ) 



(20) 



As is seen from the spectral representation, one can extract the 
quasi-particle weight from G in the limit of large imaginary 
time r at zero temperature: 



Gcr !(T (0;T — > oo) 



T=0 



tJ{q\Cqa\0)Heis. | ■ (21) 



Here, \q) t j denotes the ground state of a single hole with mo- 
mentum q in the t-J model, \0}Heis. the ground state of the 
Heisenberg model, and Et J and E^ els - the corresponding en- 
ergies; L is the linear size of the system and d its dimension. 
The last factor in (|T]) represents the quasi-particle weight, 
and we are interested, in particular, in its size dependence. 
Our treatment of the Green's function is based on a path- 
integral representation for G a ^ a i (0; r). As in the previous 
Section, we generalise our co nside rations to arbitrary spin 
S and use the representation ( Aldj ) of the Fermion opera- 
tors as a product of a Grassmann variable and an SU(2)- 
matrix. Noting that G a . a > is diagonal in the spin-indices, 
G a (0; t) ~: S a cr i S G(t), we arrive at the representation 



G(t) = 



j v[g]G UMe [g . T] U(r)g Q (p) 



I 1 



/ T>[g] e~fo dTC n\am 

(23) 

Here, G 1 is the r = r' = element of the hole Green's 
function in matrix representation, calculated for a fixed con- 
figuration of the S't/(2)-field g that describes the spins, 



G 



l-hole 



for] 



T ■ 



- Jo dT ' h l9] 



00 



(24) 



This expression results as the solution of the corresponding 
equation of motion. T is the time ordering symbol and h[g] 
is defined as the matrix (on the spatial lattice) which connects 
the Grassmann variables 77*, fj* and rj, fj in the quadratic form 
in £i_hole> c f- Eq. ([P7|). In the derivation of G(r), Eq. (j23|), 
we have restricted not only the action but also the integra- 
tion over the Grassmann variables to zero-hole and one-hole 
terms; in addition, we consider only the leading order in the 
limit T — > 0. Effectively, this implies the omission of any 
correction proportional to the fugacity of the hole, e - DJS / T ^ 
here, DJS is the creation energy of one hole in a perfect an- 
tiferromagnetic background, cf. Eq. (|J|). 
The general result for the Green's function, G(t), Eqs. ( ^3[ 
^4j), for one hole in the t-J model is valid under the conditions 
stated in connection with the derivation of the effective action, 
cf. the text after Eq. (ff^): in other words, it is valid whenever 
there is local antiferromagnetic order with a sufficiently large 
correlation length so that lattice effects play no role. 
Hole- and spin- degrees of freedom are coupled in G(t), 
because the hole Green's function G l ~ hois {g; r] depends on 
the time- and spatial fluctuations of the spin configuration, 
cf. Eqs. dl8|, fij^). A further evaluation of G(r) in the pres- 
ence of this coupling would require the exact calculation of 
G l ~ hoie {g; t] for each individual configuration of the SU(2) 
background field g x (r). This is an impracticable task. There- 
fore, to be able to proceed we have to take recourse to approx- 
imations. 

We start from the expression Eq. ( ^3| ) for G(t). In the limit 
t — ► 00, G(t) has the spectral representation (cf. Eq. ( pTj ) ) 



G(r 



00 



1 



e~ E ^ T Z q (5) 



(25) 



G(r) 



x>M Vr(rK(o) [gi(r)g r (0)} 



II 



ffdr + H«|w)} _ (22) 



G is independent of r, so that we can choose r = which 
we consider as the origin of the A-sublattice. Following the 
steps that led to the effective action in Section II, we decom- 
pose the lattice into plaquettes and integrate over the SU (2)- 
field M x which is defined within each plaquette. Next, we 
restrict the effective Lagrangian to zero-hole and one-hole 
terms, £ n \ am + £i-hole> see above. Performing then the in- 
tegration over the Grassmann variables rj r and 77*, we finally 



where Eq(S) and Zq(S) are the hole energy and the quasi- 
particle weight of the hole, respectively. 
If we neglect the fluctuations of the spin configuration com- 
pletely in Eq. (^3j), 5x(t) = const., Eqs. ( p4[ [jj|) yield 
G i-hole[ 5 = const . ;T ] = e -DJSr^ W A Eq{S) = DJS and 
Zq(S) = 1 in this case. If we would neglect merely the cou- 
pling to the spin fluctuations in G 1_hole , we would obtain a 
non-trivial result for Zq(S). In this case, the factorization 
of the hole- and the spin-part of G(r) are reminiscent of the 
spinon-holon decomposition observed in one dimension. 
The amplitudes of the spin fluctuations are controlled by the 
parameter 1/5*. In the sequel, we shall study the Green's 
function G(r — > 00), Eq. (|5j), in an expansion in pow- 
ers of this parameter, i.e., in the semiclassical expansion: 
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Eq(S) = DJS + 0(5°), Z q (5) = 1 + OiS- 1 ). The tech- 
nical details of this expansion are presented in the Appendix 
|b| We are interested in the quasi-particle weight only. There- 
fore, we disregard all corrections to the hole energy beyond 
the lowest order. This is achieved by scaling r = f/S and ne- 
glecting all terms proportional to powers of f which would re- 
sult from corrections to the lowest order of Eq(S) in Eq. ( p5| ) 
in a systematic expansion in 1 / 5. Then, our procedure yields 
an average of just the quasi-particle weight over the Brillouin 
zone: 



Z{S) = 



q 



Z 



q (5) = \-\z 



(i) 



1 



(2) 



(26) 

The fact that we cannot resolve the q-dependence of the 
quasi-particle weight within an 1/5 expansion, is easy to un- 
derstand: as it is shown in the Appendix ||, the expansion 
in 1/5 also implies an expansion in the ho ppin g amplitude 
t. Without hopping, h ab = h ba = in Eq. @, £i_ h ole be- 
comes local in the hole fields, and thus, the hole energy Eq(S) 
in Eq. ( ^ ) becomes dispersionless in leading order in the 1 / S 
expansion so that we arrive at the average. 
Our explicit results for the first two coefficients Z^ 1,2 > of the 
1/5 expansion of the quasi-particle weight in one and two i 
dimensions will be presented next. 



IV. RESULTS AND SUMMARY 

We start with the discussion of the one-dimensional case. The 
fc-sums in the results for Z^ 1,2 \ Eqs. (B7, B8) in the Ap- 



pendix are performed as integrals over the Brillouin zone 
with a cutoff at small momenta, k > fco = 2n/L. Then, Z^ 1 ' 



diverges logarithmically with the system size L: 



Z<$ = l + -=j- —In- 



(27) 



Here, we have neglected terms which remain finite in the limit 
L — > oo since they depend on the choice of the cutoff. For 
5 = 1/2, the quasi-particle weight Z of a single hole is 
known to vanish algebraically with increasing system size L, 
Z ~ L~ 2X (see Refs. |, [l0|). ThiSphphavior of Z reflects An- 



derson's orthogonality catastropheEJO: the states Cq CT |0}jj e .„. 
and |q)t,/ (see Eq. (pi|)), i.e., the bare one-hole state and the 
true one-hole eigenstate of the t-J model, are orthogonal in 
the thermodynamic limit. The orthogonality catastrophe is an 
effect of the quantum fluctuations, and therefore, it should not 
occur in the classical limit 5 — > oo. Consequently, we ex- 
pect that lims^oo X = 0. This suggests that the expansion in 
Eq. d26| ) should be reformulated as an expansion of In Z(S) 
in powers of 1 / S. Then, we get in first order for the exponent 
X the result X = (1 + U 2 / 'J 2 ) / '{4m 5) . 
At first, this looks quite satisfactory: the exact value for 
5 = 1/2 at t = is X = 3/16 (see Ref. [U)|). As t in- 
creases from zero, X also increases. This tendency is compat- 
ible with another exactly known result: at the supersymmetric 



point J/t = 2, where the t-J model is exactly solvable, one 
finds value X Sus y = 1/4; see Ref.|. 

However when we try to improve the lowest order result by 
including the next order 



we find in leading logarithmic order 



a In 



L 



(28) 



\nZ(S) 



, At 2 , 
-(1 + ^) 



— In — 

2tt5 n 4a 



— In — 

2tt5 n 4a 



(29) 

Here we have omitted a term of the order of izp{t/ J) 4 , since it 
does not change the picture qualitatively. Obviously, this ex- 
pression for \nZ(S) cannot simply be interpreted as the first 
two terms of an expansion in powers of 1/5 of the exponent 
X(S) in a power law Z(S) ~ L~ 2XI - S \ At first glance, one 
might think that the deviation from such an expansion can be 
attributed to the fact that in Eq. (p6|), the quantity Z(S) is an 
average over the Brillouin zone of the quasi-particle weights 
Z q (S) so that one cannot expect to find a power-law depen- 
dence of Z(S) on the system size L. However, for t = 0, 
i.e. for a static hole, for which Z q (S) is independent of the 
wave number q, Eq. ( p9| ) can still not be interpreted as an 
expansion in powers of 1/5. The infrared divergent terms 
in Z(S) show that instead of a straightforward expansion a 
renormalization group treatment is necessary. This is quite 
plausible, since 1/5 is the coupling constant of the non-linear 
tr-model £ n \ am , Eq. (|l6), on which our 1/5 expansion is 
based. However, in Eq. (23) the parameter 1/5 appears not 



only as the coupling constant of cr-model C n \ am , it also oc- 
curs in the hole Green's function, Eq. d24|). Therefore, the 
standard JEnormalisation group procedure for the non-linear 
(7-modelEj cannot be applied in the present case. The develop- 
ment of an appropriate procedure is beyond the scope of this 
paper. Notwithstanding these complications, Eq. ( p9| ) shows 
a behavior compatible with a vanishing quasi-particle weight 
in one dimension. 



In two dimensions, our general expressions Eqs. ( |B7| , |B8J) for 
Z\ 1 ' 2 ) yield a different picture: the k-sums converge at small 
momentum k. Thus, in contrast to the one-dimensional case, 
the average quasi-particle weight Z(S) is independent of the 
system size. This means that up to the order 1/5 2 of our large- 
5 expansion there is no sign of an orthogonality catastrophe in 
the two-dimensional model. In recent, very elaborate numer- 
ical studies of t-J model on the square latticeElEl it has been 
found that within this model, the spectral function of a single 
hole shows the signature of a coherent quasi-particle, i.e., the 
quasi-particle weight Zq(S) has been found to remain finite 
throughout the Brillouin zone. Although the semiclassical ex- 
pansion does not allow us to determine the weights Zq(S) at 
individual q points, our result for the average weight Z(S) is 
consistent with these numerical findings. 
In summary, in this paper we have presented a new approach 
to the problem of hole dynamics in the t-J model. We use 
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a generalisation of this model, which has originally been de- 
signed for spin i particles, to describe particles with arbitrary 
spin 5. This has allowed us to consider the case of large 5 
in which the semiclassical approximation is applicable. This 
approximation has the advantage that it works independent of 
the spatial dimension of the system we wish to investigate. 
Therefore, we have been able to study the single-hole dynam- 
ics in one and in two dimensions on the same footing. In 
agreement with most previous investigations of the dynami- 
cal properties of a hole in the t-J model, we find a qualitative 
difference between these properties in one and in two dimen- 
sions: while in two dimensions the results of our 1/5 expan- 
sion are compatible with the picture of the hole as a coherent 
quasi-particle in the sense of Landau's Fermi liquid theory, the 
hole appears to have a vanishing quasi-particle weight in the 
one dimensional case so that the quasi-particle picture fails in 
this case. 



APPENDIX A: COHERENT STATES 

We use the method of coherent states (for a recent review see 
Ref. |28j) ia-arder to deal with the problem of excluded double 
occupancyEl First, consider the case of spin 1/2. Then, the 
Hilbert space consists of a fermionic state with spin up cj|0), 

a fermionic state with spin down c{|0), and a hole |0) at each 
lattice site. Following Ref. [l9| we introduce (at each lattice 
site) the states 



\l>) = >c-cV] c t |Q) 



(1 - -77*77) ct |0) + 77 



with 



c f = c{ 5 TT + c{ 



.9 J [ 



(Al) 



(A2) 



Here, the parameters 77, 77* are Grassmann variables and de- 
scribe the hole, g is a 5[/(2)-matrix, g = 1 2 ; it can be 
parametrized by three Euler angles and describes the orienta- 
tion of the spin. In this Appendix, we index a unit matrix with 
its dimension. 

Next, we generalize the coherent states to the case of arbi- 
trary spin 5. Then, the Hilbert space consists of (25+1) + 25 
states at each lattice site which are the states of a spin 5, 
|5,m = -S---S) and a spin 5-i|(5-i),m = -(5- 
1 ) . . . (5_ 1 )). Now define, cf. Ref. [27|, 

\w) = (1-^77*77) R \S,S) +VR\(S-\),(S~^)) ■ (A3) 

Here, the rotation matrix R can be represented for arbitrary 
spin in the following way by three Euler angles ip, 9, and 



R 



e -i4>S z e -i9Sy g-iV-S* 



(A4) 



Instead of using the Euler angles, one can also parametrize R 
by the elements of the 5f/(2)-matrix g which is the 5 = i 



representation of R. In this paper, we use g as the fundamental 
variable. In the case 5 = ^, (A3) coincides with (Al). 



For the application of the method of coherent states, we need 
to specify the expectation values of the operators in the Hamil- 
tonian with the states \u>). Obviously, (cj||o;) = 1. The spin 
operator S acts as usual on the spin-states; with s = 5 or 



S = 5-|: 



(s,s| R 1 S R | a, a) = s n . 
In terms of g, the unit vector n reads 

n = \ Tr { <T 9<Jzg ^ ■ 
For the expectation value of the spin with \u>), we find 

(cu\S\u) = \ (25 - 77*77) n. 



(A5) 



(A6) 



(A7) 



The Fermion operators c a (we use a = T, J. = ±1) link the 
states with spin 5 and 5— |, 



c a |5, m) = 7<r,m |5-i,m - cji) 



(A8) 



Here, 71, _s = 7-1, s = and the remaining coefficients j a ,m 
are uniquely determined by demanding that the relation be- 
tween the Fermion operators and the spin is as in the case of 



5 — 5, 



1 ( t t 

2 V C T C T ~ c l c l 



S z and etc 



(A9) 



(A10) 



(All) 



(A12) 



Our convention J dg = 1 for the invariant measure of g leads 
to 



We get 

Hcjw) = V2S 77* g„\ ■ 
The resolution of unity reads in our case 

dw = l iS+l . 

The integration measure is 

dw ... = [ dg [ dr!*dii 25 e'^ 71 * 11 



dg R \s, s)(s, s\ R/ = 



1 



11, 



2s + 1 2 ^ 



(A13) 



With the aid of Eqs. ( \3A12A13), one quickly confirms that 
the dimension of the Hilbert space comes out correctly: 



Tr 



doj\uj)(uj\ 



= / d^di] 25 e~^ r > 71 (1 - 277*77) = 45+1 . (A14) 
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To complete the list of relations involving the coherent states, 
we quote the expression entering the kinetic part of the action 
(we consider now uj to depend on r) 

(w\d T \w) = (2S- V * V )(^d T g) n + r,*d T r, - ±d r (<Q*r,) 

(A15) 

(A16) 



In deriving Eq. (A.15), we used 



(a,a\tf d T R\a,a) = 2s (g^d T g) 



1 1 



APPENDIX B: 1/5 EXPANSION 

Here, we describe the simultaneous expansion of G(r), 
Eq. (p3[), in the fluctuation amplitudes and in the hopping in- 
tegral t. The antiferromagnetically ordered state, n x (r) = e z 
is given by g x (r) = 1 (cf. Eq. (^)). Then, it is convenient to 
characterize an arbitrary state by 



e i<re(x,T) e i 



e z (x,' 



(Bl) 



by a SU(2)/U(1)- 



where we parameterize the field g x 
symmetric part and a U(i)-f actor. 

It is easily seen that the latter disappears from our consid- 
erations: n x (r) is determined only by the two fields in the 
vector e = (e x , e y ) which enter the SU(2)/U (l)-part of the 
field g; the [/(l)-part, the last factor in (Bl ) drops out. Thus, 
the [/(l)-factor also drops out of £ n \ um - It still appears in 
£l-hole' Eq. (11^), through the field 



A l**( T ) = cr9 T e 2 (x,r) 



and through the off-diagonal matrix elements h^ +s x 



CJCJ 

(B2) 
and 



^x°x+5 ( c f- Eq. (19)) in which the [/(l)-terms occur as the 
first and the last factor. Looking at Eq. (|l7|) one recognises, 
however, that these U(l) factors can be absorbed into the 
Grassmann fields by the gauge transformation 



?7x(t) 



„ie z (x,T 



(B3) 

so that they disappear from the integrands in the numerator 
and in the denominator on the r. h. s. of Eq. (^3|). Further- 
more, in the integration measure e z and e separate, dg oc 
de z G? 2 esin(2e)/e, so that the integrations over e z cancel be- 
tween the numerator and the denominator of the r. h. s. of 
Eq. d23|). These considerations show that the field e z disap- 
pears completely from the expression Eq. ( |23| ) for the Green's 
function G(r). 

Now, we consider this expression in a systematic expansion 
with respect to e, i.e. we expand on the r. h. s. of Eq. ( |23[ ) the 
integration measure, the factors G 1_hole and [.gj.g ] n m me m " 
tegrand, and C n \ am in powers of e. Inspection of the quadratic 
part of £ n i(j m shows that in this expansion, each additional 
power of e beyond the second order yields an additional power 
of 1 / \/S. We make this explicit by scaling the spin wave am- 
plitude e by a factor of (45)~ 1/2 , e = e/ (45) 1/2 and the 
imaginary time by a factor of (2dJS)~ 1 , r = f/(2dJS). 



Then, the quadratic part in the new fields becomes indepen- 
dent of 5 and the spin wave dispersion reads 



/ 1 d 

= Vl-7k> with 7k = -^Y^ cos ( k sa) ■ (B4) 



"4j 



Relative to the leading order in S, G(t) = exp(— f/2), we 
want to obtain corrections up to and including the order 1/S 2 
accurately. Therefore, we retain terms up to the order e 4 in the 
factors G 1-11016 and [<?o.g ] TT m tne integrand, anharmonicities 
up to the order e 4 in £ n i CTm , and terms up to the order e 2 in 
the measure (note that after the scaling, the measure already 
contains a prefactor of 1/5 in the order e 2 ). In performing the 
integrals over the various combinations of powers of e, which 
occur after this expansion in the numerator and the denomina- 
tor of Eq. ([23l), use can be made of the linked-cluster theorem. 
G contains spatially diagonal, h aa and h bb , as well as 
spatially off-diagonal terms, h ab and h ba . The latter are of 
the order of e ~ 1/5 1 / 2 . Since they are proportional to the 
hopping integral t (cf. Eq. (jl9|)), our expansion in 1/5 be- 
comes necessarily also an expansion in t. We shall neglect 
terms of order higher than t 2 . 

The actual expansion is straight-forward but rather tedious. 
We express the result as 



G{r) = e-^ 2 (l^G«(f)-iG( 2 )(f) + 0(i)) 
For the first order coefficient, we obtain 



(B5) 



1 . . 1 _ p-^kf 



. At 2 . 
+ ^ 



(B6) 



The first term in Eq. (B6), proportional to f, belongs to the 
expansion of Eq(S) in Eq. ( p5[ ) and will consequently be ne- 
glected. In the evaluation of the momentum sum in the second 
term in Eq. (|B6|), we consider a finite lattice of M lattice points 
of the A-sublattice and exclude the point k = from the sum, 
since the corresponding mode is a Goldstone mode, the uni- 
form rotation of all spins. Then, this term has a finite limit 
as r — > 00 and defines the first order correction of the 
quasi-particle weight in Eq. (^6|) as 



dJ2 



(B7) 



The calculation of the next order coefficient, G^ 2 ^(f), pro- 
ceeds exactly in the same way, but it is much more tedious. 
We neglect polynomial terms f, f 2 , since they belong to the 
expansion of Eq(S) in Eq. (|25|), exclude the point k = from 
the momentum sums, and take the limit r — > 00. Then finally, 
collecting all terms, we get the following result for the second 
order correction in the quasi-particle weight 



1-^1 1 

A/" 2 4w k 4w k / 



.4 



k,k' 



4^ 
dJ 2 



B 



k,k' 



(B8) 
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with 

£?k,k' 



1 Wk^k' 



= -1 



4 cj k + w k / 



+ (w k + W k / 



In the preceeding expressions for A k k / and -B k , k s we have 
neglected terms of order k 2 , k' 2 , or kk' and higher. 



(B9) 



7k+k' - 7k 7k' 



Wk^k' 



1 c^kOV 
2i4 + Lj k - 



7k+k' - 7k 7k' 
WkWk' 
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